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Abstract 

Angular dependence of gap, seen in photoemission, its evolution with doping and 
temperature are interpreted on base t-t'-U Hubbard model in which a pseudogap is 
a working function for electrons removing from dielectric segments of zone bound- 
ary. Tunnel spectra in one-particle approximation display only superconducting gap 
and too large asymmetry in regions aside from optimal doping. Spectral functions 
confirm sharp well defined Fermi boundary at k x = k y direction and smoothed de- 
structed one at the T(0, 0) - M(ir, 0) - Y(ir, it) path. Both tunnel and ARPES data 
seem require spreading a model to bilayer ones. 

PACS numbers: 74.72.-h, 74.25. Jb, 74.25.-q, 71.10.Fd 
1. Introduction. 

Up to now a set of decisive experimental results have been obtained concerning elec- 
tronic structure of cuprates (see reviews [0, H, U). They include evidence of d-wave super- 
conductivity, discovery of anisotropic normal state pseudogap and "small" Fermi surface 
(FS) in underdoped (UD) materials accompanied possibly by destruction of FS 0, ||. 
Direct measurement of superconducting (SC) gap are provided by tunneling spectroscopy 
f(|- [f|. Most recent detailed photoemission (ARPES) studies reveal exotic temperature 



dependence of gap (T^, |TT| and the doping dependence of its anisotropy in superconducting 
state [Q. It was revealed the apparent discrepancy between this data and characteris- 
tics obtained from penetration depth data. It was shown that angular dependence of 
superconducting gap A(<£>) extracted from ARPES data cannot be adjusted by simple 
form ~ cos(2<y9) and the shape of gap function depends on doping. In particular, the 
ratio ( = «a/A of «a = |<iA(<£>) / d<p at if = it /A to maximum value of gap A decreases 
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essentially with doping decrease in UD materials. Finally standard shape of the FS has 
been recently revised. It seems to be large partially electron-like FS around T = (0, 0) 
instead the arc around Y(ti,tt) fll3| . 

The aim of our paper is to study which of new experiments can or cannot be explained 
in simple quasiparticle interpretation of pseudogap origin. 

Interpretation of phase diagram and pseudogap have recently |T^|-||18|| been proposed 
on base of different approaches for description of the correlation mechanism of SC pairing 
in cuprates. One of them is the spin-polaron approach in t-t'-J model |^4[ . Another is 
variational approach incorporating the correlations of the valence bond (VB) type in t-t'- 



U Hubbard model fllj^, |16| . In each case the Hubbard splitting of whole band arises due 
to long-range antiferromagnet (AF) order and fine details of upper edge of lower Hubbard 
band (LHB) determines the low energy phenomena. The hole attraction in d-channel is 
shown to be induced by spin-polaron interactions in t-t'-J model or equivalent correlated 
hopping interactions due to VB formation in t-t'-U model. According theories [^4]-[18 



the next nearest hopping t' plays important role in formation of pseudogap and a sign of 
t'/t determines its anisotropy. For systems with t'/t > a crossing of van Hove singularity 
(VHS) by the chemical potential /i corresponds to optimal doping and is accompanied by 
change in the Fermi boundary topology. The hole pockets around k = (7r/2,7r/2) with 
" small" Fermi surface are replaced by large Fermi surface in course of transition from UD 
to over-doped (OD) systems. At 5 < 5 opt the VHS occur below \i and some segments of 
magnetic Brillouin zone (MBZ) boundary near ku = (^,0) become dielectric segments. 



According [16(] the work function /i — Ek > for electron removing from this regions of 
k is just a normal state pseudogap observed in ARPES. In difference from standard VHS 
scenario |3] with united band here we deal with VHS of lower Hubbard subband. 

The method and outline of present study are following. We use the variational method 
to construct explicitly the correlated state and to test if it could display the observable 
properties. Starting from AF state with the VB correlations we explain the supercon- 
ducting (SC) and pseudo- gap anisotropy and its temperature and doping dependencies 
(Sec. 2). Then we calculate the tunnel spectra of system and discuss its distinction from 
the observed ones (Sec. 3). Manifestation of FS in one-particle spectral functions and 
A(k,u — > 0) and the FS shape are discussed in Sec. 4. 

2. The doping and temperature dependence of gap anisotropy. 
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Calculations have been done for t-t'-U-V Hubbard model by method proposed in [ 16 
Here t, t' are the nearest and next nearest hoppings; U and V are the on-site and nearest 
sites Coulomb-like interactions. The variational correlated state \I> = W(a)<& with the VB- 
type correlations (the band analogue of the Anderson's RVB states) is taken as a unitary 
transformed uncorrelated state $. Unitary operator W(a) = exp{a J2< nm >a( c na c ma — 
h.c.){n ng _ — n ms _)} , a = —a, optimizes the covalency of < nm > bonds. A problem with 
Hamiltonian H in basis of correlated states is reduced to a problem with effective 
Hamiltonian H = W^HW « H + a[H, Z] + ^-[[H,Z], Z] in basis of uncorrelated states 
$. Thus the effective problem H is solved in the mean field (MF) approximation with 
subsequent minimization of energy over variational parameter a (a < 0.22 at U/t < 9). 
Most general uncorrelated state $ of BCS-type and with double magnetic unit cell allows 
a testing the possible AF and SC pairing in system. The terms of correlated hopping 
like ~ a[/4 ff c m(J n n2: n m£ in H caused by VB formation are responsible for SC order of 
d-symmetry compatible with AF order. Band energies E k \ (eigen values of linearized H£) 
corresponds to self consistent $ and optimal a. 

The studied homogeneous solutions display the 2D AF spin order in a wide doping 
range 5 = |1 — n\ < 5 C ~ 0.3 which is larger than range of the bulk AF order (5 C ~ 0.05). 
This implies that finite radius of AF correlations is much larger than lattice constant a 
and we can use the solutions and a picture of split band as a background in discussing a 



low energy phenomena. In normal state a LHB has the following structure [16 



E k = [e + At'c x c y + ...] - y/[Udo + ...] 2 + [2t(c x + c y ) + ...] 2 (1) 

Here only main contributions from lowest harmonics with c x ^ = cosk x ( y } are explicitly 
written and d =< (—l) n S zn > is the staggered spin. In d-wave superconducting state 
the self consistent quasipartilce energies are very close to form 



E kX = ±^(E k - /i) 2 + Wl W k ~ F k {k 11 (c x -c y )w 1 + ...] (2) 

Here W k is determined by SC order parameters wi = (Z 2 — ly) < c n ^c nJtl ^ >, |/| = 1, A 3 
and F k is determined by structure of one-particle states of LHB. First harmonic gives main 
contribution to W k and its angular dependence is close to W k ~ (c x — c y ). 

The doping dependence T c (5) for Hubbard model with U/t = 8, t'/t = 0.05 and 
V/t = 0.1 is presented in Fig. la. (Note, that at V=0 the value of T™ ax increases by 
factor 1.64 |T6[-) Curve 1 in Fig. lb presents the value Wm(S) = W{k M ) of gap function 
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Wk at k = kM — (tt, 0) for the d-wave superconducting pairing calculated at temperature 
T/t = 0.002. But according (2) the quasiparticle energy and the corresponding shift 
A = Au; of the edge of photoemission distribution function at k = (ir, 0) is determined 
not only by Wk, but also by the LHB energy \Ek — n\ relative to chemical potential. 
Unlike original unsplit band, the LHB energy Ek is periodic in magnetic Brillouin zone 
(MBZ) and goes through a maximum on MBZ boundary at any path connecting the points 
T(0, 0) and Y(ir, ir). This determines the different doping dependence of ratio £ = v a/Am 
in UD and OD presented in Fig.lc. Here v& = -dA/dip at ip — 7r/4 and A M is the A at 



In OD regime each path T — Y crosses a "large" Fermi-boundary: Ek — /i = at 
k = kF and Ek — \i > at all values k at MBZ boundary. In this case a gap A(<p) in 
excitation spectrum for k changing along any path r — Y is determined only by d-wave 
SC gap Wk at hp- Therefore the gap anisotropy A(ip) coincides with A (ip) = Wk(<p) ~ 
z(<p) = 0.5(cosfc x — cos k y ) at k F . Corresponding ratio £ = v&/A M measured at k — k M 
in |T^| is equal to £ = \(dWk/ dip) /Wm ■ It is close to constant value ~ 0.7 expected for 
pure d-wave SC order without pseudogap. 

In UD systems with t'/t > a "small" Fermi boundary around the hole pockets at 
k = (7r/2, 7r/2) is formed whereas the segments of MBZ boundary near M(jr, 0) become 
dielectric segments on which Ek — \x < 0. A shift A(&m) of photoelectron distribution 



function in ARPES at k M = (n, 0) is equal to A(A; M ) = J A 2 N + W^- Here Wm is 

Wk at kM an d the value A^ = \Ek — fi\ at ku is the normal state pseudogap equal 
to work function for the electron removing from k = k]<^ in normal state. Value Ajy 
increases with doping decrease. For the same reason the anisotropy of total shift A(ip) = 



sj A 2 N (ip) +Wk((p) of photoemission edge at the dielectric segments of MBZ boundary 
differs significantly from anisotropy of Wk(<p) only. In UD region (Ajv 7^ 0) one obtains 
£ ~ 0.5[W'(7r/4)]/^A 2 N + W M , where W'(<p) = dW k /dtp at <p = vr/4. Value £ sharply 
decreases (Wk — > 0) with doping decrease in full accordance with experimental behavior 
of £(<5) [|T^]. In contrast to £(5), a quantity £w($) — v a/Wm with pure superconducting 
gap function in denominator remains almost constant (curve 2 in Fig.lc). 

Fig. Id presents the results for anisotropic shift A(<p) of photoemission edge for strongly 
UD system (5 = 0.14, T c /T c maa: = 0.55). Dependencies A(tp) are taken at k moving from 
k x = k y to M(tt, 0) along generalized Fermi boundary consisting of nonshadow part of 



k = (vr,0). 
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small Fermi boundary and dielectric segments of MBZ (see fl6f). Several features of 
calculated anisotropy A ((p) have a correspondence with experimental anisotropy features. 
1) Deflection of A(ip) from linear dependence A(<£>) ~ z((p) have really been observed 
in UD Bi-Sr-Ca-Cu-0 (BSCCO) |12] and the observed sign of next harmonics in it is 
the same as in calculated curves. 2) The region of gap A& near its nodes occurs more 
sensitive to temperature than at k ~ (7r,0). A shift A(</?) at <p ~ 7r/4 disappears with 
T — > T c whereas the pseudogap of normal state in region k ~ (tt, 0) retains nonzero 
value A M = A N = /i — E(ir, 0) > at T > T c . This feature also corresponds to 
experimentally observed different temperature sensitivity of A(</?) at different <p [|1C], p7lf| . 
Finite /c-resolution in ARPES may be a reason of more smooth observed dependence 
A(<£>) in compare with calculated one. Besides, the latter is obtained in neglecting the 
energy dispersion over k z . 3) The calculated ratio rj = 2A^ /kT c have very large values 
in underdoped case. Thus for 5 = 0.16 (T c = 0.9T c ma:r ) it is r\ — 9.5. At the same time 
for optimal and over- doping the whole gap is determined only by superconducting part 
A = W(kp) and the ratio is about rj w 5. This value is close to rj — 4.5 obtained from 
empirical BCS models with d-wave pairing. Decrease of 2A/kT c with doping increase is 
consistent with experiment 

Independent estimations of gap anisotropy near the nodes have been obtained in |~2 
from behavior of supernuid density p s at low temperature on base of relation between v a 
and penetration length at T —>■ 0. This relation dp s (T)/dT ~ d(\~ 2 )/dT ~ 1/i>a gives 
the doping dependence of v& different from the mutually consistent dependencies va(S) 
obtained from ARPES data and our calculations. Origin of distinction is unclear and 
question arises whether the temperature dependence of 1/ A 2 at T — > may be controlled 
not only by low energy Fermi excitations near the nodes of SC gap, but also by the 
low frequency spin excitations which adiabatically change the overall spectrum of Fermi 
system. Important role of spin wave fluctuations has recently been confirmed in |ZIJ . 

3. Tunneling spectra for AF+VB correlated states 

Tunneling spectra can provides a sensitive test for the obtained simple homogeneous 
solutions ott — t' — U Hubbard model which take into account the AF and VB correlations. 
Recent reliable experimental results in tunnel spectroscopy 0-0 give the base for such 
testing. 
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Consider so-called SIN (superconductor-insulator-normal metal) tunneling typical for 
STM measurements. Total tunnel current /(V) = e(Wi—Wn) is determined by probabili- 
ties Wi(u) of electron transfer from superconductor to normal metal and vice versa. These 
probabilities are found in first order in tunnel interaction V tun ~ J2a, n G{n — n )(c| lt7 c? e(T + 
h.c). Here d ea are the Fermi operators of normal metal (a tip of STM). A "tip formfactor" 
G(n — no) characterizes the extension of surface-tip interaction centered at some site tiq. 
As usually the constant density of state (DOS) for normal metal (tip of STM) and the 
energy independent formfactor are assumed. 

For state with AF+VB correlations and d-wave SC order the quasiparticle operators 
Xk\ = b\iUi\, i, X = 1, ... ,4, are eigen operators of linearized effective Hamiltonian Hi 
16|| . Here a basis set of Fermi operators is b' ki = {c' k ^, c~ , c_fc|, c_^}j, k = k — (vr,7r), 
and k G F varies inside MBZ \k x ± k y \ < it. With taking into account the structure and 
spectrum Ej^x of quasiparticles the differential tunnel current takes the form 

= W^ F { [RU ' X + kU ^'^ - e ^ + ^ + Au !x\f'(E kX + eV)} (3) 

k 

Here R = G 2 (k), R = G 2 (k) and f'(E) = df /dE is derivative of Fermi function. Deducing 
of Eq.(5) includes the averaging over the position of "tip center" no among two magnetic 
sublattices of AF state. Sign of V is chosen so that negative V corresponds to occupied 
levels in the sample spectrum. The matrix elements Uix are the analogues of the coherence 
factors in simple BCS theory. As in [^T| they are retained in explicit form in Eq.(5) since 
the normal state DOS of superconductor sharply depends on energy near VHS. Note, that 
here we deals with VHS of LHB instead of VHS of whole unsplit band as in ET . 



Figures 2,3 present the tunnel spectra (3) for nearly optimal doping and for UD and 
OD systems at temperatures lower and higher T c for one of studied variants of "tip 
formfactor" . In complete accordance with expectation at temperature lower than T c the 
SC gap (not pseudogap) displays in differential tunnel current. The distance between 
peaks asides main dip at V=0 corresponds to gap 2Asc*(<5) whose values and doping 
dependence at small T are close to dependence 2Wm{§) i n Fig. lb and is characterized by 
maximum at optimal doping. 

In OD region the calculated gap behavior 2/S.sc{&) is hi accordance with gap behavior 
measured by STM in BSCCO 0. In both cases the gap decreases with doping at 5 > 5 opt . 
In UD samples the observed tunnel gap seems to include the pseudogap which does not 
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vanish at T > T c || [![]. This result is not reproduced by calculations. Besides, at all 
doping the observed spectra retain almost symmetrical form close to calculated spectrum 
for optimal doping (Fig. 2). However a tunnel current dl/dV becomes highly asymmetric 
in UD and OD systems. At low temperatures sharp peak in function dl/dV at V < 
or V > 0) for UD and OD doping correspondingly reflects the VHS in the normal state 
DOS. The asymmetry of observed spectra remains almost invariable at slight variations 
of doping arround optimal value || [F]. The shape of DOS determines also asymmetric 
nonlinear background of tunnel current at large V (see inset in Fig. 2). However the one- 
particle approximation (3) used is unsuitable at large V. It cannot describe secondary 
wide peaks in observed tunnel spectra at large V. These peaks (as well as similar patterns 
in ARPES data) are connected by many authors with spin- wave excitations . Possible 
origin of VHS broadening or splitting removing too sharp asymmetry may consist in 
interlayer coupling of CuOi planes. We neglect this coupling. The split of VHS in bilayer 
system and position of fi between two VHS might manifest as two almost symmetric peaks 
in tunnel spectra which retains as pseudogap even at T > T c . 

Additional foundation for hypothesis on the bilayer-split have been obtained recently 
from dependencies of spectral functions A(k,u) on photon energy in ARPES data [ |13|j . 
The observed spectral features are discussed in terms of united band without huge gap 
between Hubbard subband which are implied in our calculations. One must also keep in 
mind possible inequivalence of Cu0 2 planes closest to surface and those in bulk caused 
by uncompensated fields of the charged BiO plane at cleavage in STM experiment. 

4. Manifestation of Fermi surface and influence of spin fluctuations. 

One of open question in correlated variational approach [15] was how to concert the 
ARPES data with predicted shape of nonshadow FS for OD models. For such models 
"large" FS around point r(0,0) with electron- like segments near (tt, 0) was obtained. 
This was in contradiction with widely recognized shape of FS as an arc around Y(ir, 7r) 
0. But recent new reexamination of FS shape for BSCCO shows the possibility of the 
electron-like FS crossing on line T — M(n,0). Besides the FS shape, it is important to 
discuss how various sections of FS should manifest in spectral functions. 

The second question is about a role of spin fluctuations. Important influence of low 
frequency collective spin-wave on Fermi spectrum of correlated system is wide spread 
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known (see [20| among last publication). Since the spin- wave states are the states with 
a spiral spin structure, the comparative study of the normal state DOS and the band 
structure features for AF and spiral states of t — t' — U — V Hubbard model have been 
carried out here. This study has only qualitative character. For this reason and for sake 
of simplicity we consider and compare the mean-field (MF) solutions with AF and spiral 
spin structures but without correlations of the VB type. 

Standard MF spiral state with a spirality vector Q is the MF state with one- 
electron averages < cL-c^ >= b Q exp [iQn] . At Q — > Qaf = (t 1 "? 71 ") it comes to AF 
state. Neglecting of VB correlations significantly changes the energy of system and shifts 
to larger value a critical doping at which the MF AF (or spiral) solutions transform to 



paramagnet ones [16]. For the same reason minimization of MF energy over Q at given 5 



leads to overestimated difference between optimal vector Q(S) and Qaf — (7r, 7r). In this 
connection we study and compare the DOS of spiral states at arbitrary vector Q, which 
does not minimize the MF energy at each doping. 

For MF states without VB correlations and SC pairing we define under- and over- 
doped states as that for which the VHS occurs lower or higher than the chemical potential 
ji. Then at "optimal" doping the VHS is just placed at chemical potential. For calculated 
model with parameters U/t = 8, t'/t = 0.05, V/t = 0.1, such definition gives a value 
5opt ~ 0.25 for the "optimal" doping in case of AF state. It is larger than value S op t = 0.19 
for the same system in the AF+VB state which takes into account the VB correlations. 

Fig. 4 presents the DOS of LHB for MF AF and spiral states with different directions of 
Q and \Q — Qaf\ = 0.3. The single VHS for AF state transforms to split VHS's for spiral 
states. Formation of spiral spin structure removes degeneracy of VHS at k — (±7r, ±7r) in 
the same manner as the lattice distortion removes it in usual Jahn- Teller effect || . Spiral 
states provide examples of electronic analogues of Jahn- Teller effect. 

It is instructive to compare the influence of spin structure on occupation number 
n k = \ J2a < c la c ka > m momenta space and on the one-particle spectral function 
A(k,u) at low frequency uj — > 0. 

! (4) 



A(k,u = 0) = ^EEI< fM > \ 2 eM~Ei/kT)- 



ZN^%> — " '<K[{E i -E f -uY + 1 *\ 



uj=0 



Here Ef are energies of initial and final states, Z is partitition function of system, 
and standard 5-function 5{Ei — Ef — uj) is replaced by Lorentz function with width 7. 
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Fig. 5 presents functions n& and A(k, to = 0) calculated for k changing along the counter 
y(7r,7r) — r(0, 0) — M(tv,0) — Y for under-, optimally and over- doped MF states of 
t - t' - U - V model at 7 = 0.0031 

Several features, partly well known, of obtained spectral functions must be stressed. 

Even in AF state, for which the band energy are periodic in MBZ, functions n^, A(k, u) 
do not possess such periodicity. Neglecting fine details, one can see that occupied {rik > 
0.5) or empty (n^ < 0.5) momenta states refer to regions around r or Y correspondingly, 
just in same manner as it takes place for noninteracting system. This property is well 
known and retains for spiral states also. 

At T — Y direction the most sharp fall of is seen at k x = k y < n/2 and corresponds 
to crossing of nonshadow Fermi boundary with line V — Y . Position of fall corresponds to 
the standard observed Fermi boundary in form of arc with center at Y(tt, 71"). Involving 
the spiral (instead of AF) spin order does not change the position of Fermi boundary 
at this direction. Corresponding peak in A(k, uo = 0) is most prominent at any doping 
and does not influenced by spin fluctuations, i.e. its position does not depend on the 
spirality vector Q. This is in accordance with the observed well defined Fermi surface at 
r — Y direction for any doping [|]]. The second peak of A(k, uo = 0) on the same direction 
and corresponding step in rik refers to shadow Fermi boundary. The intensity of peak is 
much less and its position significantly depends on Q. So this feature is expected to be 
smoothed by spin fluctuations. 

At any doping in region near M(tt, 0) the occupation number changes smoothly 
with k running along V — M — Y and low-frequency Fermi excitation peaks in A(k, u = 0) 
are sufficiently intensive in whole region near M at any doping, especially at "optimal" 
doping (Fig. 5). This is the same in both cases: at UD when the region near M become 
dielectric one (E(tt, 0) — // < 0) and at OD when the region near M arranges between 
main and shadow large FS's (E(tt,0) — /J > 0). Similar smoothing or destruction of FS 
in region k ~ (tt, 0) has been observed in ARPES data [||, |5]]. Large summary intensity of 
photoemission from this region of k have really been observed in |T3j for BSCCO. However 
it occurs visible only at the photon energy 22 eV, but not at 33 eV. In our picture of 
split Hubbard bands the VHS of LHB is more strong than that expected from tight- 
binding unsplit band. So additional increase of intensity at k ~ (71*, 0) is possible. But 
the singlelayer models of cuprates cannot describe the observed dependence of additional 
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photoemission on the photon energy ||13|| . As in discussions about tunnel spectra, we come 
to necessity of extending the study on bilayer models. 

In conclusion, fine structure of upper edge of LHB in t-t'-U Hubbard model and vari- 
ational correlated state with VB formation can describe some of low-energy phenomena 
in cuprates. Doping and temperature dependencies of gap anisotropy obtained from pho- 
toemission data fl2| are explained by interplay of gap and pseudogap which is working 



function for electron removing from dielectric segments of generalized FS. It is shown 
that only SC gap, but not pseudogap, is seen in tunnel spectra. At T — > differential 
tunnel current for optimal doping is similar to observed one, but it displays too sharp 
asymmetry in heavily UD (and OD) regimes with peaks at negative (or positive) voltages 
caused by normal state VHS. Model seems to need extension on bilayer models. Occu- 
pation number and spectral functions calculated in MF approximation confirm the 
sharp Fermi boundary at k x = k y direction and smoothed (destructed ) Fermi boundary 
near (-7T, 0) as it observed in ARPES @, [| . Change of FS topology at optimal doping 
and large partly electron-like FS at OD case in our model seem to receive confirmation 
in reexamined ARPES data Ul3| , but single-layer model cannot explain the dependence 



of photoemission spectra on photon energy. 
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Captions to Figures 

Fig. 1. 

a) Doping dependence of SC transition temperature T C (S) in unit [O.Olt]. b)Double 
shift of photoemission edge 2A M (S) at k — (ir,0) (curve 1) and corresponding SC gap 
2Wfc(7r, 0) (curve 2) in same units, c) The ratios £ = va/Am (curve 1) and = 
VA/Wk(n,0) (curve 2) as function of doping, d) Angular dependence of gap A(ip) for 
various temperatures in heavily UD model 5 = 0.14, T c /T™ ax = 0.55. Rest parameters 
of model are in text. 

Fig. 2. Differential tunnel current for nearly optimal doping 5 = 0.2 (T c /t = 0.0114) 
at temperatures T/t = 0.002, 0.014, 0.018. Insets: the same and the current I(V) for 
T/t = 0.002, 0.018 in large voltage regions. Parameters of model are U/t = 8, t'/t = 
0.05, V/t = 0.1. 

Fig. 3. Tunnel current dl/dV at T/t = 0.002, 0.014, 0.018 for UD and OD models 
with T c /T™ ax = 0.743 and 0.63 correspondingly. Model parameters are same as in Fig. 2 

Fig. 4. DOS of LHB for the MF AF state (curve 1) or for spiral states with Q 2 = 
(0.97r,7r) and Q 3 = (0.93tt, 0.93tt) (curves 2 and 3). 

Fig. 5. Occupation number (solid lines) and A(k,u = 0,7) with 7 = 0.003t (dotted 
lines) for MF AF state and spiral states of UD and OD models. Vector k runs along path 
Y — T — M — Y. For spiral states vectors Q are determined by values q' = 0.8n and 
q = 0.86. Dashed dotted curves on segments T — M — Y (upper grafics) correspond to 
AF state at nearly "optimal" doping 5 = 0.3. 
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